Plasma effects on resonant fusion 
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We investigate the effects of plasma interactions on resonance-enhanced fusion rates in stars. 
Starting from basic principles we derive an expression for the fusion rate that can serve as a basis 
for discussion of approximation schemes. The present state-of-the-art correction algorithms, based 
on the classical correlation function for the fusing particles and the classical energy shift for the 
resonant state, do not follow from this result, even as an approximation. The results of expanding 
in a perturbation solution for the case of a weakly coupled plasma are somewhat enlightening. But 
at this point we are at a loss as to how to do meaningful calculations in systems with even moderate 
plasma coupling strength. Examples where this can matter are: the effect of a possible low energy 
^^ C -'r^'^ C resonance on X-ray bursts from accreting neutron stars or on supernova lA simulations; 
and the calculation of the triple a rate in some of the more strongly coupled regions in which the 
process enters, such as accretion onto a neutron star. 

PACS numbers: 97.10.Cv, 26.50.-l-x 



I 

o 



> 

O 
(N 
(N 

rn 

O 



X 



1. INTRODUCTION 

Nuclear fusion rates play a central role in stellar evolu- 
tion theory. Their calculation depends primarily on lab- 
oratory data, and the open questions that persist have 
mostly to do with these measurements and their extrap- 
olations to very low energies. Another issue that presents 
itself is that of the influence of the surrounding plasma 
on the fusion rates. This is a mature subject, to say 
the least, for the case of non-resonant fusion, but there 
are some questions that still are not settled, and some of 
these pertain to the special case of resonance-enhanced 
fusion, which is the subject of the body of this paper. 

To set the context, we need to begin with a review 
of the way the plasma corrections work for the ordinary 
case in which the energy variation of the vacuum value 
of [cross-section X velocity] is mainly from the Coulomb 
wave-functions, rather than from a resonant peak. 

A lore has developed, the "basically classical" ap- 
proach [l|- [6|, which follows the following path: 

1.) calculation of an effective two particle correlation 
function (down to a small distance) in a classical simula- 
tion; 

2.) defining an effective two body potential propor- 
tional to the logarithm of this function, written in the 
form Ves = Z'^e'^/r + Vsc{r) where Vgc is referred to as 
the "screening potential" ; 

3.) in the potential Vcs calculating the tunneling 
through the barrier to determine the correlation func- 
tion at zero separation K{r = 0). This correlator then 
is assumed to embody all of the effects of the plasma on 
the fusion rate. 

Implementing the above, we consider fusing particles 
A and B, with respective number densities, n^ and rih. 
We write the rate, w, at which one specified nucleus of 
type B undergoes fusion as, 

K{0) 



and we consider all of the energy dependence of the vac- 
uum cross-section a to come from the Coulomb poten- 
tial. The factor i^o(O) is the correlator in the absence 
of the plasma, and serves to cancel the Coulomb factor 
that comes in a. We should emphasize that everything in 
the present paper pertains to cases in the thermo-nuclear 
regime, rather than in the pychnonuclear regime. 

Turning to the resonant case, in the absence of plasma 
the contribution of a narrow resonance to the fusion rate 
is, 
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where Er is the energy and F,- is what the partial width 
of the resonance in the channel of the fusing nuclei would 
be in the absence of the Coulomb force between the decay 
products. The final factor in ([2]) is the exponential fac- 
tor in the Coulomb wave-function; and /u is the reduced 
mass of the fusing particles. The parameter ^ is easily 
calculated from the development of section 2: it does not 
enter the ratios that are here addressed. 

This paper discusses approaches to incorporating 
plasma effects into the rates for the resonant case. These 
plasma effects can be important in applications. Refs. [7| 
and [8] suggest that the existence a resonant state in the 
12q_|_12(^ system in the vicinity of Er =1.5 MeV in en- 
ergy might explain features of X ray bursts from accreting 
neutron stars. Refs. [7| and [9[ discuss the screening cor- 
rections to resonant rates, motivated by this possibility. 
Following these references, we choose for our examples a 
density of p = 5 x lO^gc"^ for a pure ^'^C plasma, and a 
temperature region 4 < Tg < 9. 

For purposes of exposition we assume that the screen- 
ing potential Kic is nearly constant in the classically in- 
accessible region for the two particle system with reso- 
nant energy, and make the replacement Vsc{r) « Kic(O) 
in this region. In the above application this is a good 



approximation-'^. We write K!c(0) = 
dimensionless parameter 7 from ref. 



7T and take the logio[w/wo] 
ameter 7 from ref. [9[, 

7 = 2.293r°-25 + 1.053r - 0.5551 log(r) - 2.35 , (3) 



where T is the usual plasma coupling. Under our sam- 
ple conditions of composition and density we have F = 
.62/Tg. A first guess as to the form for the screened rate 
formula could be, 
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where the penetration factor has been altered by sub- 
tracting the screening potential at the origin, — 7T, from 
the energy in the incoming state. In normal fusion, where 
there is an integral over energy and the other energy de- 
pendent factors are the phase space and the statistical 
factor exp[i<^/T], a simple shift of integration variable by 
—^T gives the usual enhancement factor exp ji^ in rate, 
but this does not apply to the resonance case, it appears. 
In refs. [7| and [9|, the result dH) is extended by includ- 
ing a plasma-induced resonance-energy shift AE'^ given 
by making the replacement, Er — > -E^ + A^;^, and leading 
to, 



Wscr-l-shft = naCe 



-{E^+Ae^)/T 



exp 
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The question now is, "what is AErT\ Both ref. J7| and 
ref. [9| choose it to be given by £'12 — 2Ee, where Ez is 
the free energy shift of an ion with charge eZ due to in- 
teraction with the plasma. In this case, in the "basicah 
classical " approach it is a theorem that A^^ = Vsc{0) 
But in the next section we shall argue that the energies of 
the incoming ^C ions probably should not have been sub- 
tracted, in which case we would have had A^;^ ss 2Vsc{0)- 
What do we mean by "probably" ? Actually, only that 
if we were to start from what we believe to be the cor- 
rect ab initio approach, and ask how, making reasonable 
assumptions, we can get anything like the above struc- 
ture, then we would come to the latter conclusion. When 
we look more deeply into the problem, however, it looks 
more and more as though nothing about this whole gen- 
eral approach can be trusted. 

For now, just in order to show the importance of the 
issues, we compare the two suggestions made above for 
the parameter range specified above. In fig.l we compare 
the ratios of screened to unscreened rates for our two 



^ It is easy to go beyond this simplification, using the mechanics 
presented in ref. 
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FIG. 1: The logarithms of the ratios of the screened rates, 
wi and W2 to unscreened rates wo from ([5]). The rate W2 has 
been calculated taking Ae^ as just the resonance energy shift, 
while the rate wi is the result if, in calculating Ae^ we also 
subtract the energy shifts of the incoming nuclei. 



different assumptions as to the value of A^^, , plotted over 
a range of temperatures. 

Looking at fig. 1, we can stipulate that the subject of 
this paper can be important, even in a region of moder- 
ate values of the plasma coupling. But we are not here 
addressing data, nor the likelihood of existence of the 
possible resonance. We are only examining the underly- 
ing theory, and we must characterize ([5]) as an ansatz, 
rather than as a derived result. This is in contrast to 
what we would argue vis-a-vis the case of most applica- 
tions of ^ to non-resonant fusion, but even the latter is 
subject to the following caveat: eq. (P), giving the rate 
correction in terms of K, is true only in the case that the 
energy release Q in the fusion is much greater than the 
Gamow energy. Eg [10| UJ. For the case of resonant 
fusion the analogue to the value of Q is, in effect, not 
only not large compared to Eq, it is negative. It is our 
opinion that we can only begin serious consideration of 
plasma effects in the resonant case only after developing 
a correct formal framework for the discussion. This is 
taken up in the next section. 



2. FORMAL FRAMEWORK 

We quote some general relations that are derived in 
the appendix. In these derivations for the fusion reac- 
tion, A+B— >R, we take the resonance, R, as a particle, 
with a point coupling to the fusing nuclei. Described 
in terms of field operators tpa{j^) and '0b(r) for the fus- 
ing particles, and '0r(r) for the resonance, a zero-range 
fusion interaction is. 



Hfus — . 



dr tpl ipa ipb + H.C. , 



(6) 



where calculation of the decay rate for R^-A+B gives g^ 



in terms of T 
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Since the above model gives back exactly the Brcit- 
Wigner formula when we sum the chain of repeated inter- 
actions to find the scattering amplitude, we judge it to be 
adequate for our purposes, as long as the resonance width 
in the fusion channel is small. The total width will not 
enter the explicit results, but for the whole application 
we assume that the decay is mostly into other channels. 
This simulation of the resonance should be good as long 
as the barrier penetration calculation doesn't care about 
the difference between r — Q and the nuclear radius. 

For purposes of exposition we think of singled-out re- 
acting particles, A and B. The "plasma" is composed of 
all other particles. The rate will then have a factor V^"^ 
where V is the volume; when we restore translational in- 
variance and have finite densities of the species A and B 
this factor is to be replaced by the product of number 
densities naUh to get rate per volume, or just by Ua to 
get rate per nucleus B as in our convention stated above. 
We make this replacement ah initio in what follows. It 
will also eliminate superfluous complication from what 
follows to take both the nucleus B (described by V'b) and 
the resonance R to have infinite mass, and to position 
them at the origin. The corrections, expressed as a mul- 
tiplicative factor, will be very nearly the same as in the 
equal mass case. 

With this simplification, in the appendix we derive ex- 
pressions (|A.8|) and (jA.lip . which together give the reso- 
nant fusion rate, 
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where /? = l/T and Zp is the partition function for the 
plasma. We have used the notation, PI, to stand for all of 
the coordinates of the plasma particles collectively. The 
symbol B within the first bracket indicates that nucleus 
B is present (at r = 0) in the space of states in which this 
calculation is to be carried out. In the second bracket we 
have only the plasma and the single nuclear state R, the 
latter fixed at the origin. 

Now we ask: "What is the relation between (|8|) and the 
ansatz ([5]) ?" We begin by ignoring the coupling of the 
resonance state R to the plasma. For the final bracket in 
([8]) we have simply. 



i?, PI' e^'*" R, PI 



^-iE,.t 



5{V\ - PI') , 
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where the delta function is the product of delta functions 
for the individual plasma particles. Then the sum over 
plasma states in ([8|) gives the dependence exp(— i^E^O 
times K{Q,t), where K is the correlator as used in ([T|), 
but now extended to finite time by the replacement 
/3 — > /3 — ii and in the form of a trace over the plasma 
coordinates PI. In the "basically classical" approach it is 
just. 



K{0,t) 



1^ 
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where the VPp is the wave-function at the origin in the 
screened potential. For the case considered in the intro- 
duction in which the screening part changes little within 
the classically inaccessible region we obtain. 



l*p(0) 
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Putting (dni), ® into dH) and using ([TT]) we obtain 
(HI), later to be modified by the plasma induced energy 
shift. Next we ask, "Does the governing formula ^ tell 
us how to do this modification"? Note that once we 
put in interactions with the plasma the final factor in 
dS]), {R,Pl'\...\R,Pl) will not be a delta function in the 
plasma coordinates, nor diagonal in these coordinates. 
Thus we cannot eliminate the sum over states PI'. More- 
over the first bracket is no longer a trace and is then no 
longer given terms of the correlation function. 

However, in the spirit of "basically classical", we can 
say "We determined the effective potential between ions 
A and B as the logarithm of the classical correlation func- 
tion, in a calculation in which traced over the plasma 
states. So we now imagine that when A and B are close 
to each other the plasma is in a single configuration, PIq, 
a configuration that produces an electrical potential 4>sc 
at the origin. The sum of the interaction energies of 
these two particles with this potential, less the Coulomb 
energies of the individual widely separated ions, gives the 
screening potential energy V^c(O)." 

With this state PIq replacing both PI and PI' in ([8]), 
and eliminating the sum over PI', the second bracket then 
becomes. 



i?, Plo 
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which leads back to ([5]). 

With this pseudo-derivation in hand, we feel more con- 
fident in declaring our belief that if any approach that 
leads to the form (O makes sense, for strongly coupled 
systems, then A^^ must be simply the energy shift of the 
resonance R in the plasma, as given by e^{Z\ -\- Z2)(j>sc, 
with no subtractions for the energy shifts of the isolated 



A and B states in the same plasma. Indeed, the object 
estimated in (|12p has no way of knowing how the total 
charge Zi + Z2 was split between the incoming nuclei. 
The above conclusion is at variance with the assump- 
tions of refs. [7| and [9|, and we have already illustrated 
in fig. 1 the very big difference that the choice makes in 
the rate formula. 

That said, we believe that it is much more likely that 
the basically classical approach, which appears to have 
served us well in ordinary fusion reactions, when the en- 
ergy release Q is much greater than the Gamow energy, 
simply cannot be used to get meaningful results in this 
resonance fusion problem, with any choice of A^^ . 

We best can explain the essential difference between 
the resonance case and the normal fusion case by recall- 
ing some of the features of the latter. Again, strictly for 
economy of exposition, in the reaction A -I- B — >■ C -I- D, 
we take the two ions B and D to have infinite mass. As 
in the resonance case we take a zero range fusion inter- 
action. Instead of a resonance energy parameter we now 
have an ordinary energy release parameter, Q. The ana- 
logue to (j5]) is now, as given in ref. jll|, 



.g^Zp^n^/r^die^Q* 
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Now we can see how in the ordinary fusion case, when 
Q >> Eg the "basically classical" approach may suffice. 
In the A-(-B — >-C-|-D fusion process of (|13l) when the en- 
ergy release Q is large compared with the Gamow energy, 
it is a reasonable approximation to set t = in the expo- 
nent within the first bracket () (pertaining to the initial 
A + B system). If one further omits all plasma coupling 
of the fusion products to the plasma in calculation of 
the second bracket (pertaining to the C-l-D system), and 
keeps just the kinetic energy term, Hq, then the time in- 
tegral just sets the energy of the fusion products to the 
value Q. Since we have taken the operator in the second 
bracket to be independent of the plasma coordinates, we 
can integrate 6{P\ — PI') over the intermediate plasma 
configurations PI', and the first (), integrated over the 
space PI, becomes the trace defining the static correla- 
tor, bringing us back to the standard result (HJ. 

But for the resonance case we find no similar argument 
that we can make to defend p^ or to lead back to ([S]). 



3. PERTURBATIVE RESULTS 

As an alternative to the machinations of the last sec- 
tion, we can, beginning from ([8]), turn the perturbation 
crank to generate results for the weak coupling case. 
Though the results would of course not be applicable to 



the strongly coupled systems discussed earlier, we could 
hope that they would cleanly resolve the question of what 
to use for the resonance energy shift in ([5|), assuming 
for the moment that the calculation leads to the general 
structure of ([5]). But it will not, as we shall see. Nonethe- 
less, in the appendix we do work out these corrections in 
detail. For one thing, they comprise the only certain 
result as to the plasma corrections in the resonant case. 
For another, they highlight problems that were swept un- 
der the rug in the discussion that led to ([T2|) . and that 
raise doubts about any efforts to make a connection to 
the "basically classical" lore, as used in our speculations 
above. 

First we give the result. The expansion is in powers 
of the couplings of the plasma particles to each other 
and to the fusing particles. The e's in the Coulomb 
wave-functions of the fusing particles are not part of the 
expansion, of course. Because of the long range of the 
Coulomb potential, the leading corrections are of order 
e^KD, where kd is the Debye momentum. 

After a long calculation, given in the appendix, involv- 
ing the cancellation of many terms of this order, we find 
simply. 



w = [1 + e'^ZiZ2KDl3]wQ , 



(14) 



which is just the original Salpeter correction, as used in 
standard weak-coupling applications, most particularly 
in the sun since that is the only venue in which relatively 
small corrections are of much interest. Is this what we 
did so much work to learn, that the effect of the plasma 
is still just in the initial state equilibrium contact prob- 
ability? That is, whether the final system is resonant or 
not, we have the same modification factor? Indeed if we 
look back at the ansatz , (|H) and choose A^,. = — 7T, 
as in the case in which the screening potential does not 
change appreciably within the classically inaccessible re- 
gion, we recover the equivalent result, namely just the 
same screening factor exp(7) that would modify a non- 
resonant fusion ^. 

However, in the perturbative calculation in the ap- 
pendix we see that the simple result depends completely 
on including the terms which, in the notation of ([5]), the 
intermediate plasma states PI' are different from the ini- 
tial plasma states in the trace, PI. Described in terms of 
fundamental processes we can say that the incoming nu- 
clei scatter a plasma particle and then this plasma parti- 
cle is put back into its initial state through an interaction 
with the final resonance; or in an equivalent description 
in terms of scattering amplitudes rather than correla- 
tion functions, that we must calculate interference terms 
between amplitudes in which a plasma particle scatters 



^ This was pointed out in ref. Q. 



from a fusing nucleus and amplitudes in which it scatters 
from the resonance. 

If we left out terms like this and calculated only the 
separate contributions from the interaction of the plasma 
with the initial state and from the interaction with the 
resonance, then the sum of these terms would not be 
given by (|14l) . Depending on the parameters, it can be 
an order of magnitude greater than shown in (jl4[) . Yet 
going back to the strong coupling case, combining the 
two effects, classical screening and energy shift, is exactly 
what was prescribed. So it appears to us that it is not 
meaningful to calculate this way, and that, for example, 
neither of the "energy shift" models that we used in the 
comparison of fig. 1 are viable. 

For comparison, if we did the same kind of pertur- 
bation calculation for the normal fusion case we would 
again obtain the Salpeter result with additional terms 
that cancel when we add them up. In contrast to the 
resonance case, however, these individual terms are all re- 
duced by at least one order of the factors {ZxZie^nu)IQ 
or £^Gamow/Q, as long as these ratios are small. Thus 
in this case the perturbation analysis gives us no cause 
to question the qualitative argument for the "basically 
classical" approach given at the end of the last section. 
In ref. Ill[ | we have argued the need for rethinking the 



plasma corrections for cases in which Q is so small that 
the condition is not fulfilled. 



4. DISCUSSION 

Our strongest conclusion in this paper, which lacks 
definitive computational results for the systems of most 
interest, is that one should begin with a complete and 
correct formulation. It appears to us that the literature 
on this topic has depended on using physical intuition to 
graft together some pieces of lore that have served well in 
other contexts (but with some limitations that are usu- 
ally unacknowledged). We do not believe that this pro- 
vides a coherent approach starting from basic physics. 

We list the two new results of this paper that have 
been proven, then turn to speculation: 

1). The formula for reaction rate, ^. This was stated 
for the case in which one of the fusing particles and the 
final resonance both have very large mass, so that they 
each would stay at the origin. But it provides an ade- 
quate testing ground for the general questions that we 
raised in relation, e.g., to ^^C -|-^^ C fusion. '^ 



2). The summation of the leading correction terms 
in the perturbative development. The outcome argued 
against our favorite guess as to how to do the energy 
shift calculation, but the details did show the necessity 
of following the plasma in a way in which one does not in 
the "basically classical" approach. We characterize the 
procedure of the latter as, "Do one big classical simula- 
tion, for the effective potential, and thereafter do your 
fusion calculation with the individu al p lasma states out 
of the picture." There is evidence [12] that this proce- 
dure is adequate for determining the contact probability 
in the initial state, for an ordinary fusion reaction, as 
long as -EGamow << Q- But when Q is small, or we are 
dealing instead with resonance fusion, we believe that it 
is insupportable. 

If we had to choose between the two alternatives of- 
fered in the first section, we would follow the discussion 
just before and after ([T2|) and use the energy shift of the 
resonance particle itself for A^^ , with no subtractions for 
the energy shifts of the incoming particles. This gives the 
large relative suppressions that we see in the \og[w2/'Wo] 
curve in fig. 1. 

Although in the above opinion we have reverted to the 
kind of guesswork that we deplored in the introduction, 
we do add that the perturbative results do not argue as 
strongly as may have appeared for using the other op- 
tion for Ae^, as adopted by refs.[7] and [9,]. First, sup- 
pose that we separate out the the normal Salpeter fac- 
tor, (1 -I- e'^Z'^KoP), from all other terms. Now the sum 
of these other terms is convergent even if we calculate 
them with unscreened Coulomb interactions, i e., the in- 
frared divergences in the individual terms have cancelled 
against each other. 

This behavior is similar in form to the cancellation 
of (photon) infrared divergences in the inclusive cross- 
sections for charged particle scattering, where individual 
contributions, like photon bremsstrahlung, are divergent 
in an order by order expansion. The point we make is 
that these cancellations do not carry forward in any fash- 
ion to the non-infrared-divergent parts. Returning back 
to our perturbative calculation, we believe that the can- 
cellations of the infrared terms that restored the Salpeter 
form do not extend in any sense to the non-infrared di- 
vergent terms (which begin with terms of order e^ but 
also include higher corrections of order F'^). Therefore, 
for applications in which the coupling is even of moder- 
ate strength, the only thing that we really learn from the 



^ The corresponding formula for the case of general masses is 
easy to derive, but more complex. In the bracket for the fus- 
ing particles, the two particles that meet at in the left hand 
state, ( I, meet again in the right hand state, | ), but now at po- 
sition r. Correspondingly, in the bracket for the fusion product 
R we go from position r to position 0. The whole expression has 



an additional integral over r; all of this serves to transfer total 
momentum from the initial to the final state. Center of mass mo- 
tion is no longer totally trivial when we are in a plasma and Q 
is insufficiently large. The derivation of the above follows easily 
from the methods used in the appendix. For the ordinary fusion 
case the analogue is shown in detail in ref. [lCl( | , with a somewhat 
different representation of the operators involved, however. 



elaborate calculations in the appendix is that it is wrong 
to disregard the changes in the plasma, i. e. the sum 
over the states |Pr) in ^. 

We ask how one might imagine doing a real calculation, 
assuming unlimited computing resources. The guiding 
example could be the calculations of Militzer and Pol- 



lack [1^ or Ogata [13|, who do real quantum path in- 



tegral evaluations for the static correlator at small dis- 
tances, without the assumptions of the "mostly classi- 
cal" approach. In ref jl^j the authors write e~^^ = 
jg-^ff/iVjAf g^j^^ choose N sufficiently large so that each 
factor (Pl'|e~^'^/^|P1"|), effectively at high temperature, 
can be expanded in perturbation theory. The product 
of the TV factors, involving N integrals over intermedi- 
ate plasma configurations, is very expensive to compute, 
needless to say, but because of the negative real expo- 
nentials the calculation is feasible. In contrast, the result 
(fT3)) . with its factor e^'-^^**^'^ and extra time integral, 
cannot be calculated directly by such a procedure, be- 

* Our refs. Q and [9|| are cited in |15| . but their results not used 
in the analysis. 



cause of the oscillations. Clearly some technical progress 
is greatly needed. 

An important fusion chain that requires resonance en- 
hancement is the triple a process. There are sug gest ions 



for treating the associated plasma effects in ref. [IJ]. In 
red giant cores where the process makes its first appear- 
ance, the plasma coupling is too weak to make very much 
difference. But when, e.g., it takes place in surface lay- 
ers of accreting neutron stars, all of the caveats of the 
present paper apply. 

We also note recent considerations of the role of a pos- 
sible ^2C-|-^2C resonance on the ignition of a type la 
supernovae 15| [16|. Though screening issues were not 



addressed in either of these works, the plasma coupling 
is strong enough in this application that they need to be. 
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Appendix 

For completeness we begin by recapitulating some development from ref [lOJ . We calculate the time rate of change 
of the species A induced by the fusion interaction ^. Since the system is translationally invariant we can choose 
to evaluate this time derivative at at point r = and choose time to be zero as well. Directly from the Heisenberg 
equations and the commutation rules, writing the transform of © as. 



Hins{t) 



dr [h{r,t) + h{r,ty 



where h = gipaipbipl '^^ obtain the rate of change of Ua, the density of species A, 



The notation {...) p indicates the thermal average in the medium, such that for an operator, O, we have 

{0)p = Zp'Tr \0 expi-/3[H + Hf^, - fiaNa - ^ibNb]) 



(A.l) 



(A.2) 



(A.3) 



where Zp is the partition function, and where H is the Hamitonian in the absence of the fusion term. We wish to 
calculate the rate to lowest non-vanishing order; i.e. to second order in h. Thus we now must consider the linear 
response of the average of the operator [h{0, 0) — h\0, 0)], as it appears in (jA.2[) . to the perturbation Hfus, giving an 
expression for the rate, 



-I I dt /(dr)(K(0,0),i?fus(t)] 



dt / dr([h{O,0)-h\o,0)],[h{r,t) + h\r,t)] 



(A.4) 



where now the thermal average in the medium is to be calculated under the action of H alone. We note that, 

[[h{0, 0), h{r, t)]) ^([h^O, 0),h\r, t)]) = . (A.5) 



Using in addition the space-time translational invariance of the medium and the antisymmetry of the commutator, 

(^[hHo,0),h{r,t)]y=-(^[h{O,0),h\~r,~t)]) , (A.6) 

we can write the rate as 

w = - f dt f dr/[h{O,0),h''{r,t)]y. (A.7) 

When we take the medium to contain no nuclei of type R, so that there is no reverse reaction, we can omit the first 
term in the commutator in ()A.7p . Then inserting ^ we obtain. 



w 



g' f dt /'dre-'^-*(v4(r,t)V4(i-,i)V'r(r,t)V4(O,0)Va(O,0)^b(O,0)) =5' / dt y{t), (A 



where, 



y{t) = Zp^ / rfrTr[Ae-(^-")^Vl(r,0)^l(r,0)V'.(r,0)e-"^Vj(O,0)Va(O,0)7/;b(O,0)], 



(A.9) 



and where A = exp[/3^/^iiVi]. 

The relation (jA.8|) is the fundamental formal result required to determine the fusion rate. We have taken the 
Schrodinger and Heisenberg representations to coincide at i = 0, and at this time the ip operators can be expanded 
in annihilation operators for single particle states. 

The calculation that follows will be much simpler if we specialize to the case in which one of the fusing particles, 
say B, and the resonance, R, are much heavier than ion A; furthermore the plasma corrections to the order that we 
calculate, as a multiplicative factor, do not depend at all on the mass configuration. Therefore we simplify, ab initio , 
by taking the masses of B and R to be infinite, and situating a single B at r ~ 0. We replace operators ipb{^, t), '0r(r, t) 
by operators b{t), and r{t) that respectively annihilate the indicated particles situated at the origin. Note that even 
though ions B and R no longer carry kinetic energy, these Heisenberg operators are still time dependent by virtue of 
the Coulomb interactions of the ions to which they refer. The result (|A.9[) is now replaced by , 



y(t) = ZpV-iTr[Ae-('3-**)^V'l(O,0)6t(0)r(0)e-"^rt(0)^a(O, 0)6(0) 



(A.IO) 



The rate that we calculate is now the rate of production of the resonance from the single ion B situated at r == 0. 
When we recast the expression in a form in which we consider only one A particle, the factor V~^ in (JA.IOP gets 
replaced by a factor of Ua. Written a bit more explicitly to exhibit the plasma coordinates in an intermediate state 
(JA.IOP reads. 



y{t) 



Ua X 



E 
pi.pi' 



ra = 0,B,Pl 



^-{P-it)H 



ra = 0,B, Pl')(i?, PI' 



itH 



i?,Pl 



(A.ll) 



which leads directly to ^ . 

To calculate the leading terms in a perturbation expansion we write the total Hamiltonian (not including Hf^s) 
as H = Hq + Hj , where Hq includes all kinetic energies as well as the mutual Coulomb interactions of the plasma 
particles with each other and the Coulomb interaction between the distinguished nuclei A and B. Then Hj simply 
contains the Coulomb couplings between the nuclei, A , B, C and the plasma particles. Designating the total electrical 
potential of the plasma particles at position r as (j){r) we have. 



Hi = dr eanair)<l){r) + [etbH + err'' r]<j){0) 



(A.12) 



Of course, if we were attempting to go beyond the weak coupling regime, we would have to face the fact that the 
potentials ip in (IA.12P involve all of the coordinates of the individual plasma particles and that a numerical simulation, 
as will be required for a definitive answer in the strongly coupled case, will involve the explicit taking into account of 
the paths of all of these particles. In the discussion section we describe the difficulties in principle of carrying out such 
calculations. But in the present paper we calculate to second order in Hj, and in addition single out only the terms of 
order e^ in which the superficial perturbative order e^ has been lowered through the long range of the Coulomb force. 



For this purpose all we need to know about the plasma is the part of the field-field correlation function in momentum 
space that is singular for small k as kd goes to zero, namely [10|, 



0(k,t')0(-k',i") 



IR 



47r(5(k-k')/3"^K|, 



(fc2 + KJy)k-^ 



We introduce the interaction picture through the identities, 



-HAfi-^t)=e-''o(0-^'^n+{-t^,t) 



-,-iH,4 



where, 



n'-+\-iP,t) = exp 



i I dt'Hiit') 



f7_(0,i)e-'^°\ 



i / dt'H'j{t' 





rj(-)(0,i) = exp 



time-ordered and anti-time ordered, respectively, where. 
Substituting in (|A.9|) we obtain. 



(A.13) 

(A.14) 
(A.15) 

(A.16) 

(A.17) 



In (|A.17p we have also used the cyclic property of the trace to move the operators ipa{O,0)b{0) to the front; in this 
form the needed diagonal matrix elements in the trace are between states with no B or R present. We now expand 
the ft factors in powers of the coupling, retaining only terms of order e^, (or e^, Caeb, etc.). Schematically, each will 
have powers of e coming from the thermal expectation of a product of two ion-electric-potential (f){r) operators; these 
are determined by the Hamiltonian Hq. If this potential-potential correlator is itself expanded in powers of e, the 
expansion begins with terms of order e^, so that the rate corrections superficially would be of order e** but the small 
k singularity of (jA.13[) reduces the order to e^. In the present paper we pursue only these e^ terms; the neglected 
terms will go as power e"* and higher. 
Expanding, we obtain, 



and 






f7(")(0,t) = 1 + i / dtiHiih 



dtiHi{ti 



-ip 



i/3 



dt2Hi{t2), 



dhHiih) / dt2Hl{t2). 



(A.18) 



(A.19) 



Our task is to put (JA.ISP and (|A.19p into the rate equation (|A.17|) and select out the terms that are of order e^kd- 
We define ?/(t) = yi (^ + j/2 (0 + J/3 (0 where j/i , 2/2 , 2/3 come respectively from the second order piece of J7'+-'(—i/3,i); the 
product of the first order parts of J7'^+^(— i/3,i) and flfl~ {0,t); and the second order part of flf^ (0,i). The operator 
under the trace now factors into a part with a 0(k, /:')</)(— k',i") and a part that depends only on the operators in 
the space of the nuclei A,B,R. To the order that we need to calculate, the partition function also factorizes into a 
product in each space. Therfore we can directly use (jA.13[) to evaluate the plasma part of the trace. Since the electric 
potential of the plasma, </> occurs in Hj in the term / d'kna(k)(j){k) the k's in the in (jA.lSp in the plasma factor still 
link to the k's in the fusing-particle factor. However since the leading (order e^) behavior comes from very small k, 
we can evaluate at fc = in the fusing particle sector and use. 



f dkdk'U{k,t')(t>{~^k',t")\ = {2TTf(3~'^KD. 
Furthermore we can make the replacement, for example. 



(A.20) 



Hr 



eaB 



iHnt 



dpdk 



ipap-(-k<A(k)e 



-iH„t 



dp 



(2 



— a^ap+k / dk0(k,t), 



(A.21) 



since in the limit of small k, Hq commutes with J dp aj^Op+k. Then we see that in the infrared limit there is effectively 
no ti,t2 dependence of the integrands in lA.18[) and (|A.19[) : performing these integrals gives the respective factors, 



32 J.2 



and 



/3^ i^ 
f2{t) = {t + tl3)t, 



hit) 



(A.22) 



(A.23) 



With the t dependent factors separated from the matrix elements of operator products that we still have to evaluate 
we now write, 



yi{t) = (3>iDfiit)Z~^V-^(^0i„0. 



b'rr' 



}|Ofc,0,) 



(A.24) 



where 



Ob, Or) stands for the state in the B,R part of the space that contains nothing; Aa = exp(/3 yUaA^a); the trace 

(A.25) 



acts in the A+plasma space; and, 

W{t) - 2-3/27r-2M3/2/3K^y-ie^"'^ / dE Ve e-^^^^'^'' 1*^(0 
In similar fashion we obtain 



y2{t)^l3^iDf2{t)v-'z-'(Ob,o, 



~(l3~it)E^ 



XrCrT'C 



TraAa{[Epap*p(0), 

(ea + eb)ej2{t)W{t) , 



5a Z]q«qaq + eb&''& 



(A.26) 



and, 



y3W = /3^D/2(0V^"'^a"'(06,0,|TraA„{Epap*p(O)e 

• e^ [^^^]' [ Ep' ap'V (0)*] foVt } |0b, 0,) 
= e^j3{t)W{t). 



-{fi-it)E^ 



xr ( 



(A.27) 



In the above we did not impose charge conservation e^ = &a + &b, in order to best exhibit which interactions 
were contributing to which terms in the answer. Now eliminating e^ and summing the above three pieces, using 
h{t) + h{t) + hit) = /3V2, we obtain simply. 



y{t) = 2'^'^n'^Ml'^PnDV-\ea + e,)^e^-0 / dE ^ e-'^^-'^'^^ a^^aE *b(0) 



Substituting in ()A.8|) . w — g'^ J dty{t)e-K.p{—iErt), performing the time integration, then using, 

g^ - 2-i/2^-iAf-3/2p^^i/2 ^ 

and, 

e'^'^" =n, (2^)3/2 (M„T) -3/2, 



(A.28) 
(A.29) 
(A.30) 
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and replacing the factor V ^ by n;,, we obtain for the change of the fusion rate due to the plasma interactions 

Sw = -(Ca + ebfliKDWo , 

where 

/ 27r \3/2 _p IT ( 2 

wo^nanfcl— — 1 Vr-e ''' exp(-7reZ 
This gives, after the zero'th order term is added, 



2 .72 



(A.31) 



(A.32) 



= "'oil + ^I^HDi^a + ebY] ■ 



(A.33) 



Before drawing conclusions from (|A.33|) we note that the factor Ua in the zero'th order term emerged computationally 
from, 



i<, = (yz)-i^(e-^(^-^- 



-''»^")at„ 



p"P/ 



_ ^/SA'a 






-Pp^/2M 



(A.34) 



But if we set out to calculate Ua in our basic formalism we would have had begun instead with, 
n. = (VZ)-' j;(e-«''-''-~-l<>;op) = j;(e-''»"at<,pn(+l(-,:;3,0)) 



('4^-'^*'7|p' 



-fip^ J2M 



(A.35) 



and similarly for n\,. That is to say, the plasma interactions have changed the relation between number density and 
chemical potential in just such a way that when the relation (|A.33p is rewritten in terms of the actual density, the 
terms of order e^ (and e^) are removed [lO|, leaving us with the final answer ^, 



= u;o[l + eae(,/3K£,], 



(A.36) 



the result (fT4|) quoted in text. 
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